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Abstract 



A study of the possible vacuum phases in a strongly interacting two flavor light quark system is presented. 

Four possible phases are found with some of their properties presented. The possible existence of a 
localized diquark condensed phases inside a nuclcon and excited hadrons is investigated by showing the 
potential of solving three selected current "puzzles" in experimental obervations if such a possibility is 
taken into account. Base on these results, the potential of producing these phases in heavy ion collisions 
is speculated. 
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Chapter 1 

A General Introduction 



The strong interaction vacuum is believed to undergo certain kind of phase transition in high tempera- 
ture and density environments. Such a behef at high temperature is substantiated by both lattice QCD 
computations and the statistical bootstrap argument Q of the early days when the fundamental 
particles of strong interaction are considered to be baryons, mesons and corresponding resonances. The 
destination phases of the transition from the present day large scale strong interaction vacuum phase, in 
which the temperature and density can be regarded as zero, and the nature of these transitions remain, 
however, poorly understood. The best known new phase of QCD is the one called quark-gluon plasma 
(QGP), which was found on lattice QCD simulations 0. Other phases, such as the chiral symmetry 
restoration phase [^, was also shown to happen together with the QGP in lattice QCD calculations. 

To understand these questions is a worthwhile endeavor since, on the one hand, the behaviors of 
QCD are not known well enough to have all its possible phases listed with their properties predicted 
and, on the other hand, the experimental probes of the hadronic matter in the laboratories and the 
analysis of resulting data have not provided us with sufficiently clear pictures for the nature of the 
vacuum phases produced. For example, it is still not known for certain whether or not the QGP has 
been created in heavy ion collisions. However, an understanding of the above mentioned problems is of 
fundamental importance since the knowledge gained in this area can help us to understand the evolution 
of the early universe, the properties of the stellar matters in heavy stars, and eventually to tell, e.g., 
where matter come from, why there appears to be an arrow for the time, etc.. It thus constitutes one 
of the frontiers of our knowledge that awaits to be explored. 

QCD is a highly nonlinear theory that is difficult to solve, especially at low and intermediate energies. 
Lattice simulations are still limited by the power of computers, which allows the computation only on 
small lattice (size) systems. In addition, fermions are not easy to handle on lattices. The approach 
adopted by this study for the theoretical explorations is to use simplified fermionic models for the strong 
interaction which inherit most of the basic symmetries of QCD at long distances. These models can 
be considered to have large overlap with the ones that are derivable by integrating the gluon degrees 
of freedom (in the path integration sense) out from the original QCD Lagrangian density. Despite 
such an expectation may in fact not easily fulfilled in reality, we still consider this undertaken valuable 
since it can provide sharp predictions that can be compared to our empirical knowledge due to the 
relative simplicity of the model approaches. In addition, the simplicity of the model approach can let us 
investigate large set of possibilities (for the fermionic sector) that is otherwise not possible. In contrast 
to the theoretical exploration efforts, the model dependency of the study related to realistic physical 
systems, in which the theoretical possibilities is searched for, is reduced as much as possible by using 
model independent methods like, e.g., symmetry constraints, sum rules, etc.. In doing so, the conclusion 
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derived can be more reliable. 



On the observational sides, the information about a nucleon in various reactions, the multiplicity 
production in e~^e~ annihilation and many others can be used to determine whether the theoretical 
possibilities discovered in model studies actually exist inside static systems like a nucleon and in a time 
dependent system like in the e~^e~ annihilation system. Albeit a negative result of the above mentioned 
investigations does not necessary mean that these possibilities are forbidden in other processes like 
heavy ion collisions, which are the main topic of this workshop, the early universe, etc., a positive result 
from these "domestic" experiments does increase the probabilities and are thus worth of studying. 
In addition, the structure of a nucleon, the hadronization mechanism in the e^e~ annihilation are 
interesting problems on their own rights. 

The report consists of two parts. In the first part, which is based on Refs. ^ and work undertaken 
Ref. 1^], two massless fermionic models for the strong interaction are introduced with their vacuum 
phase structure determined. The massless fermions are identified as up (u) and down (d) quarks. Four 
kinds of phases are found: 1) the bare vacuum 2) the so called a-phase, in which the original chiral 
SU{2)l X SU{2)r symmetry is spontaneously broken down to a flavor SU{2)v symmetry, that is widely 
studied in literature 3) the /3-phase, in which the original chiral SU{2)i x S'C/(2)r symmetry is also 
spontaneously broken down to a flavor SU{2)y symmetry; it spontaneously breaks the U{\) symmetry 
corresponding to baryon number conservation induced by a condensation of diquarks 4) the w-phase, in 
which the original chiral SU{2)l x SU{2)fi symmetry is unbroken. The U{1) symmetry corresponding 
to baryon number conservation is also spontaneously broken, like in the /3-phase. The second part 
consists of a study of three different issues. The first one is related to an examination of the PCAC 
relationship for a nucleon, which is based upon Ref. 0. Certain inconsistency that is in favor of the 
existence of a localized /3-phase inside a nucleon is revealed. The second one provides a gauge invariant 
modification of the Gerasimov-Drell-Hearn (GDH) sum rule required by current experimental data. It is 
based upon Ref. |8|. If the deviation of the experimental data from the GDH sum rule is proven genuine, 
a spontaneous broken down of the U{1) symmetry correspond to baryon number conservation inside a 
nucleon is shown to be a necessity. One of the most natural realization of such a symmetry breaking is 
by having /?- or w- phase inside a nucleon. The third one is based upon the observed baryon-antibaryon 
rapidity difference correlation in a high energy e~^e~ annihilation which is discussed by Prof. Xie's group 
Q in Shandong University in this workshop. Experimental evidence exists that baryon number 
and antibaryon number is not locally produced in the observed jets. Such a violation of locality can be 
shown to be a direct consequence of the fact that the string that fragments in to hadrons contains either 
the /?- or to- phase. In the last chapter of this report, a summary is provided and some speculations 
concerning the possibility of producing the (3- or/and uj- phases is given. 
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Part I 

Theoretical Exploration 
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Chapter 2 

Preliminaries 



2.1 Introduction 



The basic participants of the strong interaction we know of today are color carrying quarks (antiquarks) 
and gluons with basic interaction between them mediated by gluons according to the QCD Lagrangian 
density, which is a non-Abehan gauge theory with asymptotic freedom at short distances or high energies. 
At low and intermediate energies, the interaction between color sources (including quarks and gluons) 
becomes strong enough to render a perturbative analysis useless. For a light quark system, strong 
interaction can cause quark-antiquark pair production from the bare vacuum. With the increase of 
the interaction strength, it is expected that a macroscopic number of quark-antiquark pairs can be 
produced resulting in effects that survive the thermal limit to lead to various phase transitions in the 
vacuum. 



At a formal level, the gluon degrees of freedom can be integrated out to obtain a pure quark-quark 
interaction effective action. Such an effective quark-quark interaction action are expected to be very 
complicated if can be practically computed at all and difficult to analyze to extracting useful information 
without certain physical picture being built up using other indirect approaches. One of these approaches 
is to build models for the quark-quark interaction part that inherits the basic symmetries of the original 
Lagrangian density. It is expected that the major physics emerged from these models has large overlap 
with the one implied by the QCD Lagrangian. 

Before list the symmetries considered, let's make the first simplification. The behaviors of certain 
massless fermion system are expected to represent the light quark system, which is restricted to the up 
(u) and down (d) quark sector in the flavor space here, well. Thus, for the simplicity of the discussion, 
we shall assume that the mass of the u and d quarks is zero. QCD Lagrangian density with massless 
u and d quarks has, besides the ones listed in the following, a SU{2)l x SU{2)fi chiral symmetry. So 
the symmetries that these models should have in addition are 1) invariance under the Lorentz group 2) 
global U{1) symmetry corresponding to the conservation of baryon number 3) local U{l)em symmetry 
corresponding to the electromagnetism 4) global color SU{3)c symmetry. 

One of the well studied model with the above properties is the Nambu Jona-Lasinio (NJL) model 



1 11 1 used at the light quark level (some include the strange quark as well). 
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2.2 Auxiliary field method using a simple model as an example 



The fermion-fermion interaction terms are assumed to be of contact 4-fermion form in the first part 
of the report. One of the best methods of treating these nonhnear 4-fermion interaction models is to 
introduce auxiliary fields |12, 13|. The subtleties of introducing auxiliary fields related to the Fierzing T 
and Crossing C operations are oiscussed in more detail in Ref. . It shall not be discussed here. Suffice 
it will be to use the following heuristic steps to show how auxiliary fields are introduced to handle these 
non-linear terms. 



For simplicity, let's consider the following fermion-fermion interaction model 

v2 



^ {i^ — m) + Gq (ipip)' 



(2.2.1) 



where m is the mass of the fermion, ip is a 4-component Dirac spinor fermion field and Go is the coupling 
constant. 

In the path integration formalism, the generating functional W[r],fi] for the connected Green func- 
tions for the fermions is 



exp{iW[7],fi]) = N DliplDlipjexp 



(2.2.2) 



where rj, rj are external Grassmanian fields, is a constant that are of no physical effects; it is so chosen 
that I^[0,0] is zero. 

Since the fermion fields are not bilinear in the arguement of the exponential integrand, the functional 
integration over the fermion fields can not be easily computed. A step forward can be achieved if the 
following mathematical transformation is used, namely, 



exp {iW[r], r]]) 



N' I i?[V']D[V']exp 
D [a] exp 



i d x{aa + bil^ipY 



(2.2.3) 



where a is the auxiliary field introduced and a, b are arbitrary constant to be determined in the 
following. The second functional integration over a gives an ij and f] independent constant contribution 
since it is a complete square in the arguement of the exponential been integrated; it causes no physical 
consequences since its effects can be absorbed in to the normalization constant N' . After writing the 
generating functional in this form, progr ess ca n be made by a proper choice of the constants a and b. 
Let a = 1/2^0^ and b = ^/G^, then Eq. ^71^ is 



exp(iTy[r?,r?]) = N' D\(t] / exp 



i i d X [C + rjip + ipr]) 



(2.2.4) 



with 



ij) (ip — a — m) 



4Gn 



-a 



(2.2.5) 



Wh en the space-time dimension is taken to be 1 + 1, this is the half bosonized Gross-Neveu model 
|13| |. It can be seen that the 4-fermion term in the original Lagrangian density are absent in the new 
one; this allows the integration over the fermion fields be carried out. The price that one should pay 
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is to introduce additional field,, namely the auxiliary field a, to be functionally integrated. Eq. 2.2.4 
becomes 



(2.2.6) 



where the effective action Seff for the auxiliary field a after integrate the fermion fields can be expressed 
as 



SeffW] 



-iSpLn7'^ {ip — a — m) 



1 



4Go 



(2.2.7) 



with "Sp" denoting the functional trace of the corresponding operator. Formally, Seff can be expressed 
in terms of the eigenvalues A of operator D = ^^(if) — a — m), namely. 



SeffW] 



1 



(2.2.^ 



where X[a] is the eigenvalue of D and the summation is over all of the eigenvalues considered. 

In case of a independent of space-time, which is the case for the vacuum state of the system, Eq. 



2.2.5 can be simplified further to 

Seff{<T) = 



-Wa 



df^p 
(2^ 



In 



— (o" + mY 
p2 _ 



4Go 



(2.2.9) 



where V4 — > oo is the space-time volume of the system. The vacuum a value of the system can be 
determined by minimizing the effective potential Veff{cr) = —Seff{cr)/V4. 



2.3 8-component Dirac spinor 

The best representation for the Dirac spinor for the path integral formalism developed in this study is 
the 8-component "real" one |Q, ^, § as oppose to the 4-component one. An 8-component "real" Dirac 
spinor can be written as 

where ipi and ■02 are 4-component Dirac spinors. They are related to each other in the following way 

, f C-ipT One flavor o 1 1 \ 

^2 = [Cir.^pf Two flavor ^^'^-^^^ 

where C = —C~^ is the charge conjugation operator (in the 4-component representation of the Dirac 
spinor) and T2 is the second Pauli matrices acting on the flavor space of the fermion field ^. 

For later discussions, it proves useful to introduce three Pauli matrices Oi,2,3 acting on the two 



4-component Dirac spinor ■01,2- With Oj, Eq. 2.3.11 can be more compactly written as 

t = -i/^n, (2.3.12) 
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with the O matrix defined as 



n 



-c-^p-^ 

Cp 



(2.3.13) 



where p = p ^ = 1 for one flavor case and p = iT2 = —p ^ for two flavor case. 



Using the 8-component the Lagrangian density corresponding to Eq. 2.2.5 takes to following 
form 



(2.3.14) 



The effective action for a is then 



i 1 
Seff M = - -SpLn7° (ip-a-m)- — 



d'^xa'^ + const 



(2.3.15) 



due to Eq. |2.3.12| . Here ^^const" is chosen such that 5ejj[0] = 0. 
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Chapter 3 



Models 



3.1 Classification of the 4— fermion interaction terms 

The light quark system in reahty contains fermionic u and d quarks with three colors. Therefore the 
fermion fields that should be used (in the 4-component form) is 4 x 2 x 3 = 24 component. For a modeling 
of the quark system, a Dirac spinor ipfc with f = u,d labeling the flavor and c = 1,2,3 labeling the 
color has to be used. For the compactness, the flavor and color indices of ip will be suppressed in the 
following as long as no confusion is thought to occur. 

If the mass of the light quarks is assumed to be zero, the full Lagrangian density of the system is 
written as 

C = i^ip^P + dnt. (3.1.1) 

The 4-fermion interaction terms can be generally classified into two categories in the quark-antiquark 
channel 

^ I singlet g J ^ ^'""^^ 

— A'nl + A'nt' (3.1.2) 

where generates quark-antiquark scattering in color singlet channel and c[^J^ generates quark- 
antiquark scattering in color octet channel. 

For vC^IIj, the well known two flavor chiral symmetric Nambu Jona-Lasinio (NJL) ||ll| interaction is 
chosen, namely 

4nl = GoU^Pi^r + i^P^7'r^pf]. (3.1.3) 



The color octet is written in the quark-quark (antiquark-antiquark) channel form for our 

purposes, namely 

^int = } triplet } sextet + ^ ^) 
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= (3.1.4) 

The color sextet term is repulsive in the one gluon exchange case. Due to the non-existence of colored 
baryon containing three quarks in nature, it is assumed to be generally true. So we restrict ourselves to 
the attractive color triplet two quark interaction terms. The attractive color triplet quark bilinear terms 
can be classified according to their transformation properties under Lorentz and chiral SU{2)l x SU (2)/j 

(3) 

groups. In general, if only terms without derivative in fermion fields are considered, has the following 
form 

= lY^GrY^cM^Mn^), (3.1.5) 

r ah 

with r^, matrices in Dirac, flavor and color spaces generating representation "r" and satisfying the 
antisymmetrization condition 



(3.1.6) 



for the fermionic systems. 



Operator ipV^'ij^ belongs to an irreducible representation "r" of chiral, Lorentz and color groups and 

tpT^ip belongs to the conjugate representation. Coefficients C^^ render the summation J2ab ■ ■ ■ invariant 
under Lorentz, chiral SU{2) l x SU {2)ji and color SU{3)c groups. {Gr} is a set of independent 4-fermion 
coupling constants characterizing the color triplet quark-quark interactions. The tilded fermion field 

operators ip and tp are defined as 

i> = ^lJ^{-iT2)C-\ (3.1.7) 
^ = CiraV'^. (3.1.8) 



The transformation properties of a list of bilinear products of two fermion fields in color triplet and 
sextet representations are given in Table 1, where all of the possible ones without any derivative in 
fermion fields are included. 



3.2 Model 1 and its vacuum phase diagram 



The 4-fermion interaction terms are classified in the above section. The representation "2" of table |3.1| 
is chosen to form the first model interaction term. In terms of the 8-component Dirac spinor ^ and 
after half bosonization p, g] , it takes the following form 



1 



XcX", 



'if 



(3.2.9) 



where a, vr, Xc and x*^ are auxiliary fields, (x^)^ = —Xc Go and G3' are coupling constants of the model. 
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Table 3.1: Transformation properties of various bilinear fermion operators under the action of Lorentz, 
chiral SU{2)l x SU{2)r and color SU{^)c groups. 



Representation 



1 
2 



5 
6 
7 



Operators Lorentz Group SU{2)l x SU{2)r SUjSjc 



( jjTlp \ 



Pseudoscalar 
Scalar 

Pseudoscalar 
Scalar 

Axial Vector 
Vector 

Vector 
Axial Vector 

Tensoi — 

Tensor — 
Tensor + 



1 dim 
1 dim 

6 dim 

4 dim 

4 dim 
1 dim 
6 dim 



3 
3 



3 
3 



3 
3 



If the auxiliary fields do not depend on the space-time, the effective potential can be computed in 
the momentum space, it has the following form 



Veff = 4i 



d'^p 



In 



+ X 

p2 



p2 



1 2 1 2 



4Go 



(3.2.10) 



where = 



-XcX 



A numerical evaluation in Euclidean momentum space shows that the minima of Vgf f (a, x) is located 
on either the a axis (% = 0) or the x ^^is (c = 0). Fe//(cr, 0) and Fe//(0, x) are found to be 



2 1 
veff{a,0) = 3/(^) + 



^A2^ 167raoA2' 



Veff{0,X) 



1 



X 



A2 ^ 167ra3/ A2 



(3.2.11) 
(3.2.12) 



where A is the Euclidean p^^ space cutoff that defines the model, the dimensionless effective potential is 
defined by V^jj = A^Ue//, ao = GqA"^ /Att and ay = GyA'^/Sn with 



-X + In ( 1 + ^ ) - ln(l + x) 



(3.2.13) 
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The value of Vejf at the minima of Eqs. |3.2.11 and 3.2.12 determine the vacuum of the system in the 
one loop Hartree-Fock approximation for the fermions. It is easy to r epr esent the phase structure of 
the model by showing it in the aQ-ay plane, which is given by Fig. 3A. Three kinds of phases for 
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Figure 3.1: The phase boundaries between the a-phase, 7-phase and the 0-phase. The chiral symmetry 
is unbroken in both the 7-phase and the 0-phase. The a-phase breaks the chiral symmetry spontaneously 
down to a flavor symmetry. 

the vacuum are possible. The first phase, which is called the 0-phase, is the bare vacuum. The second 
phase, or the a-phase, has non-vanishing average value of the chiral SU{2)l x SU{2)ji symmetry is 
spontaneously broken down to a SU{2)v flavor symmetric one in this phase. The third phase, defined as 
the w-phase, has non-vanishing diquark and anti-diquark condensation characterized by a non-vanishing 
X^; chiral symmetry is unbroken in this phase. 

The phase transition across the boundary between the 0- and the a- phases (ag = 7r/12) and the 
one between the 0- and the u- phases {a^, = vr/S) are of second order. The phase transition between 
the a- and the uj- phases (ao > 7r/12 and ay > tt/S) is of first order. The Meissner effects for the 
electromagnetic field are present in the w-phase. This is discussed in more detail in Ref. for model 
II in the following. We shall relegate such a discussion for this model to other work. 



3.3 Model 2 and its vacuum phase diagram 



The second mod el interaction Lagrangian density chosen are constructed from the ones in representation 
"4" of table |3.1| . In terms of the 8-component Dirac spinor \I' and after half bosonization ^, its 
Lagrangian density is 



1 



a 



f7"03 + O 



(+) 



6^7VA + <^^ 



13 



4Gr 



(3.3.14) 



where 4>f,c, (/"^ with )c 



(5^ with (Jt )c = are auxiUary fields introduced. 



The effective potential Veff can be computed. A numerical evaluation of it in the Euclidean mo- 
mentum space shows that the absolute minimum of Veff{a'^, cjP') is located either at cr^ 7^ and (fP' = 
which is the a-phase, or at cr^ = and (fP' 7^ 0, which is called the /3-phase, in the spontaneous 
symmetry breaking phases. The phase diagram is obtained by minimizing 14// with respect to cfP by 
assuming cr^ = in the /3-phase and with respect to o"^ by assuming 0^ = in the a-phase. The result 
is represented in Fig. 3.2. When = or (/>^ = explicit expressions for the effective potential can be 



as 
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Figure 3.2: The phase boundaries between the a-phase, /?-phase and the 0-phase. The chiral symmetry 
is unbroken in the 0-phase. The a-phase and /3-phase break the chiral symmetry spontaneously down 
to a flavor symmetry. 



derived. They are 

K//(^^o) 



and 



^ 








i 


47r 


11 


(-!- 

Vao 


-I) 





A2~ 



2-K 



a 



r2 i A4 27r 




1/1 4\ 1 

4 Va3 ~ A2 ^ 6^A4 

4a3A2 27r I + ''a2 + 3 



i.2 < A2 
^2 > A2 



(3.3.15) 



(3.3.16) 
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with the reduced dimensionless couphng constants ag ^-nd a^, defined as 



ao = (3.3.17) 

47r 

03 = . (3.3.18) 



The value of o"^ that minimize T4//(cj^,0) satisfies 



2 



A2\ / ^ \ 



It can not be solved explicitly. The value of (fP' that minimize T4jj(0, c/)^) is 

The general form of Eqs. ( p.3.15| )-( |3.3.20 ) are obtained by the replacement 

^ cr^ + vr^ (3.3.21) 
0' ^ <^^c0^' + 4c • (3.3.22) 

following the symmetry properties of The transition is of second order across bound arie s between 
0-phase and the a-phase (ag = 7r/12) as well as 0-phase and /3-phase (as = vr/4) in Fig. |3.2| it is first 
order phase transition across the boundary between a- and /?- phase (ao > 7r/12 and > 7r/4). 

The chiral SU{2)l x SU{2)ji symmetry is spontaneously broken down to a SU{2)v fiavor symmetry 
both in the a-phase and in the /3-phase. 

Some of the other properties of the /?-phase are discussed in more detail in Refs. They are not 
reproduced here. 



3.4 Spontaneous separation of baryon number in the /3- and u- phases 

The vacuum of the models introduced are studied based on the assumption that the lowest energy state 
of the system (vacuum) contains vanishing baryon number density. This assumption can be phrased in a 
different way, namely, that the baryon number and the antibaryon number in the vacuum of the system 
cancel locally leading to net baryon number density zero at each space-time point. This assumption is 
not apparent a priori since it is not an independent one. Whether or not it is true depends, as it is 
shown in this study, upon the interaction of the system. Can the vacuum state of an interacting system 
contains net baryon density locally? The answer is yes 

To investigate this question, a Lagrangian density C differ from the original one Eq. ^.1.1| by an 
additional f^aj% term, with /i" a statistical gauge field H and the baryon number current density, 
can be used [Q]. In the 8-component representation for the Dirac spinor, it takes the following form 

C = ^^{ip + 03^)^ + CM- (3.4.23) 
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In the equilibrium statistical mechanics formulated in a path integration Language, such an additional 
term is the one needed in a grand-canonical assemble. For the (3- and uj- phases, the question of what's 
the configuration for /i° in the lowest energy state of the system, which is the one corresponding to 
the vacuum by definition, can be studied using this Lagrangian density in the conventional formalism 
treated in the Euclidean space The resulting effective potential for the w-phase as a function of 
fi = VTi^/A with a definite ay (and therefore x^) is shown in Fig. 3^. It can be seen that the values for 
fj, that minimize the effective potential are nonzero; this is somewhat counterintuitive. For the /?-phase, 
the result is similar Q. It is shown in Ref. Q that in the a-phase, the lowest energy configuration for 
/i is located at fj, = 0. 
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Figure 3.3: The effective potential v^ff = V^ff/A'^ as a function of x = /u/A. 



A nonvanishing fi implies non-vanishing local baryon number density in the vacuum of the system. 
It means that baryon rich and antibaryon rich regions can be generated spontaneously in the (3- and 
Lo- phases of the systems considered. A nonvanishing /i also spontaneously breaks the CP invariance 
of the original system. These properties of the oj- and /?- phases, in which the U{1) phase symmetry 
corresponding to the baryon number conservation is spontaneously broken down, have the potential of 
solving 0] the old problem in the idea of the matter-antimatter symmetric universe [24| based on 
standard big-bang cosmology. 



The above mentioned results are actually not difficulty to perceive. In the a-phase, correlated 
quark-antiquark pairs condense in the vacuum, they can not be spontaneously separated. In the /3- and 
w-phases, correlated quark-quark pairs and antiquark-antiquark pairs condense; it is conceivable that 
they prefer to separate locally while keeping the global net baryon number unchanged. This qualitative 
picture can actually be more rigorously substantiated by using field theoretical method. The details are 
given in Ref. 0]. 



^For a consistenf discussion of the a-, (3- and lj- phases, certain new computation method has to be introduced. They 
are given in Ref. |6|. 
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Part II 



Traces of P- and lu- Phases in Nuclear 

Systems 
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Four possible phases of the hght quark system are discovered in Part I of this report. Since the models 
introduced are not derived from the QCD Lagrangian density, their relevancy can be established either 
by deducing them from the fundamental theory (QCD) or by checking whether or not the consequences 
of these models are observable in nuclear systems, like, the static case of inside a nucleon, the dynamical 
processes involving hadronic reactions, etc. 

Deducing the quark-quark interaction terms are not the goal of this research; it is certainly an 
nontrivial and interesting topic to be studied in the future. The relevancy of the results of the first part 
of the report are assessed by looking for the possible existence of the phases found there in realistic 
nuclear systems. Two kinds of systems are considered: 1) inside a nucleon and 2) in the high energy 
e^e~ annihilation processes. 

The large scale vacuum phase for the strong interaction vacuum is in the a-phase. There are large 
body of empirical evidences supporting such a notion. The reason behind the search for the /?- or uj- 
phases inside a nucleon and in the hadronization processes of the e~^e~ reaction is inspired by the facts 
that sufficient high baryon density can cause the strong interaction vacuum to flip into f3- or to- phase 
0, ^. There might be sufficiently high baryon density inside a nucleon and in the hadronization process 
of the e^e~ reaction. 

The results derived in part I are for a uniform system with its spatial and temporal extension going 
to infinity. When these results are used in case of a nucleon which is of finite extension in space or 
in case of an e~^e~ annihilation reaction which is finite both in time duraction and spatial extension, 
fluctuations effects has to be taken into account. Since what is seeking in this work is the null effects, 
which means that the effects are absent without the existence of the /3- or lv- phase, we belief that the 
question of existence of these phases are not affected by the effects of the fluctuations and a mean field 
picture is still a good starting point for systematical improvement of the result for technical reasons. 
This will be elaborated in more detail in the following discussions, which do not depend on mean field 
pictures. 
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Chapter 4 

Extended PCAC Relation for a Nucleon 



4.1 A QCD chiral Ward identity and PCAC relation 



The QCD Lagrangian density 

L = -\TrG^'^G^u + ^iip-mo)i,, (4.1.1) 

o 

with G^'^ the gluon field strength tensor, ijj the quark field, has a global chiral SU{2)i x SU{2)r 
symmetry if mg is set to zero. There is a QCD Ward identity for the divergence of 

df^A^ = 2moiji-f^T^tP (4.1.2) 

in the mo 7^ case at the "classical" level. Here r'* (a=l,2,3) is one of the Pauli matrices in the isospin 
space. Unlike its U{1)l x C/(l)i? chiral symmetry, which is anomalous, the SU{2)l x SU{2)r chiral 
symmetry of QCD is anomalies free. Therefore Eq. [4.1.2| continues to be valid when the system gets 
quantized, which turns it into an operator equation, it has, however, no predictive power in this form 
since the matrix elements of the right hand side (r.h.s.) of it between physical hadronic states is not 
immediately known. 



The PCAC relation states that 

d^AI = mlUcPl, (4.1.3) 

While QCD chiral Ward identity follows directly from the QCD Lagrangian, the PCAC relationship 

is an empirically law, which agrees with data well when q'^ ^ ~ 0. It implies the following weak 
equality 

C ~ ijij^T^ip (4.1.4) 
for their matrix elements in the momentum transfer regions considered. 



How does one go from the QCD chiral Ward idnetity Eq. ^.l.l to the PCAC relation Eq. ^.1.6? 
How good is it when q^ is allowed to go away from the pion mass shell? Or, how much difference does 
one expect between d'^A'^ and fT^rn^cj).,^ when their matrix elements between single nucleon states are 
taken? These questions are studied in this work. 
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An necessary condition for it to be satified is to has its matrix elements between physical hadronic 
states satisfy the same equation. The PCAC relation in the mesonic sector, namely its matrix elements 
between meson states, is known to be satisfied. Is it satisfied in the baryonic sector? A closer look at 
it is necessary in order to have an answer to it. 



4.2 PCAC relation for nucleon states 

The matrix elements of d^A^^ — fj^mlcj)-,^ between single nucleon states are parameterized as 

{p' \d^^Al{Q) - Uml<t>^m p) = mlCiq^)U{p')n'T'^Uip) (4.2.5) 

with C(g2) a measure of the error of the PCAC relationship. In terms of various nucleon invariant form 
factors, it can be written as 

mgA + {q"^ -ml)gp/2 + ml [g^NN fn - mgA] = ml{q^ - ml)C{q^). (4.2.6) 

Here m is the mass of a nucleon. gA, gnNN, fir and gp are nucleon axial vector current form factor, 
the pion-nucleon coupling constant, pion decay constant and the nucleon pseudoscalar form factor 
respectively. They are functions of both q^ and m^. 

Let's define two functions 

A{q^,ml) = 2mgA + {q^ -ml)gp (4.2.7) 
B{q^,ml) = gnNNfn-mgA, (4.2.8) 

where the ml dependence is written explicitly. It follows from Eq. 4.2.6| that 



A{q\ml) = 2mlCiq^ml), (4.2.9) 
B{q^ml) = -mlC{q\ ml), (4.2.10) 

by noting that A{(f' ,ml), B(q'^,ml) and C(g^, m^) are slow varying functions of so that the coefficient 
of the explicit q'^ dependent term in Eq. [4.2.6 should vanish (for a more detailed analysis of this 
assumption, see Ref. M). 



4.3 Chiral symmetry and two basics relations between nucleon form 
factors 

It can be noted that in the chiral symmetry limit ml 

A{q^,0) = lim [2mgA + {q^ - ml)gp] = 0, (4.3.11) 

which represents the conservation of the axial vector current in that limit. It is also expected that 

lim A{q'^,ml) = 0, (4.3.12) 
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which means that the effects of ttt-q in the QCD Lagrangian can be neglected when the momentum 
transfer >> m^. In the region ~ m^, A{q'^,m'^) = m'^C {q"^ , ml) , which is small as it will be 
shown in the following. So, our first basic assumption is 



A{q\ml) 



Aiq^O)=0, 



which leads, in the realistic case of rriQ ^ 0, to the the following equation 

mgA{q^) + ]^{q^ -ml)gp{q^) = m^C(g^m^) 



0. 



(4.3.13) 



(4.3.14) 



The correction to it is of order 0(m^/M^) ~ 1% {Ma is the lightest meson next to pion in the pionic 
channel). This equation is supported by recent experimental measurements |14, ^] within < —q^ < 0.2 
GeV'^. The agreement of the above relationship with experiments is quite good. It is important further 
investigation of the relationship between gA{q'^) and gp{q'^) can be carried out. 



Here, we shall assume Eq. 4.3.14 to be true (within an error of order m^/M^ ~ mo/M ~ 1 — 2/ 
Then 



gTTNN{q^)fTT{q^) - mgA{q^ 



-mlC{q^,ml 



0. 



(4.3.15) 



The correction to it is also of order 1%. It is indeed the case on the pion mass shell q^ = where the 
value of gAirn^), g-wNNi'm'^) and f-Kim-n) can be deduced from experimental data. Data suggests that 
mlC{ml,ml) -1-2%. 



4.4 The determination of valid region 



Let's evaluate the q"^ dependence of C{q'^,m'^) so that the range of validity of Eqs. 4.3.14 and 4.3.15 
can be assessed by using an once subtracted sum rule, namely 



C{q') 



C{ml) + 



J2 _ ^2 ^oo 



vr 



ImC{s) 



sth {s-q^){s-miy 



(4.4.16) 



where the dependence of C(g'^,m^) is suppressed for simpicity. Since at small Ig^l, the details of 
ImC{s) is not important, we can use a step function to make an estimate, namely, ImC{s) « a6(s — Sth), 
with m^a ~ 1 — 2%. In this simplified case 



vr \sth- mi 



(4.4.17) 



where the dependence of C(g'^,m^) is suppressed for simpicity. Since at small the details of 
ImC{s) is not important, we can use a step function to make an estimate, namely, ImC{s) « a6(s — Sth), 
with m^a ~ 1 — 2%. In this simplified case 



a . 



C{q') = C{mi) + -ln 



Sth 



vr 



Sth - 



(4.4.18) 



In order for the corrections to our two basic equations Eqs. 4.3.14 and 4.3. 15| to increase by a nother 
1%, |C(g^) — C(m^)| has to increase by 100%, which means that the region of validity of Eqs. 4.3.14 
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and |4. 3.151 is < Q.lsth- The range of validity of these two equations in the negative region is much 
smaller than -0.7sth- 



The next step is to estimate the value of sth- The lowest value of sth in the pionic channel is below 
9m^, which corresponds to an anomalous threshold for the three pion state. However, the effective Sth 
correspond to that of the pn threshold, which is considerablly larger than 9m^. This is a consequence 
of the fact that the underlying dynamics of QCD is chiral invariant except for a small mass term. 
From this fact the dynamical (operator) equation or QCD chiral Ward identity follows. This dynamical 
equation ensures that the operator ^ij^r^Os"^ can only excite a longitudinal vector excitation since it 
is proportional to the divergence of an axial vector operator field. Therefore the state in which the three 
pions are all in a s-state is dynamically forbiden. The allowed state which dominates the dispersion 
relation is the one where two of the three pions have a relative angular momentum of 1, which is itself 
dominated by the p excitation strength. Therefore Sth ~ {Tnp + m-n-)'^ and the range of q in which Eqs. 
4.3.14 and 4.3. 15| are valid is 30 to 35 times larger than « 0.02GeV^. 



4.5 Test of PCAC relation for a nucleon 



Eq. 1.3.15 gives a specific relation between 5a(q' ) and gnNN{Q )fn{Q )■ We study whether or not it is 



consistent with the phenomenology next. 

The dependence of gAi'l'^) in the space-like region is of a dipole form jl^, namely, 

9a{<i') = T-^^^, (4.5.19) 

with Ma ~ 1 GeV. We shall use Ma = 1.0 GeV in the following. 



The q'^ dependence of g-wNNil^) is known less well than that of gA{q^)- A monopole form for it, 
which can be parameterized as 

A 2 2 

g.NN{q^) = g.NN ,^'""2 (4-5.20) 

is agreed upon in the literature; the value for A,rAfAf varies. It is found to be greater than 1.2 GeV in 
nucleon-nucleon (NN) scattering and deuteron property studies. It is in contradiction to the expectations 
of many chiral nucleon models for the nucleon. Lattice QCD evaluation also indicates a smaller one, 
namely, A^^^N ~ 800 MeV. On the phenomenological side, Goldberger-Treiman discrepancy study [|17|, 
ppTT^ v.s. pmr~^ coupling constant difference [18|, high energy pp scattering [|T9| and charge exchange 



reaction [|^, etc, support a value of A^attv close to 800 MeV. By introducing a second "pion" tt' 
with mass 1.3 GeV, A-j^nn can be chosen to be around 800 MeV without spoiling the fit to the NN 
scattering phase shifts and deuteron properties [21|. This picture was later justified by a microscopic 
computation in Ref. p^]. 

The q^ dependence of f-n-{q'^) is little known from direct experimental observations. It can be ex- 
tracted from the following time ordered correlator 

qMq')S^'' = -^(g2-m2)|dW«V)(0|^C(^)4(0)|0)(-)- (4-5.21) 



The pion is a composite particle in QCD. Eq. 4.5.21 can nevertheless be constructed using the 



following procedure. First, consider a three point function with two quark fields (in the pionic channel) 
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and one axial vector current operator. It contains a pion pole at = m^. Second, solve, e.g., Bethe- 
Selpeter or whatever suitable equation(s) to obtain the independent vertex function (or the wave 
function) for a pion. Third, project out the pion contribution to the three point function at arbitrary 
by using certain orthogonality relationF] between the pion vertex functions and the other parts of the 
quark-antiquark scattering amplitude (the four point function) in the pionic channel. Therefore the 
correlator in the above equation can in principle be computed from the QCD Lagrangian. The result 
is also expected to be unique. It is however hard in practice to obtain a reasonable result since QCD 
has not been solved. We need to resort to more controlable methods that connect to experimental data 
and is accurate enough. 

At low g^, it is beneficial to express it in terms of an once subtracted dispersion relation, namely, 

/.tf) - fAO .l^iy^^J^M^^. (4.5.22, 

This is because the resulting f{q'^) does not depend on detailed shape of /m/^(g^) but only some low 
moments of it when q^ is sufficiently small. The lightest physical state connects to the axial vector 
current operator with the quantum number of pion is the pvr two particle state, the value of Sth is 
chosen to be {mp + m-,^)'^, where nip = 770 MeV. At s' = 4m^, which correspond to the lowest invariant 
mass of a NN system, another branch cut for f-wiq^) develops. We shall include pvr state only since 
NN state contributions to the above equation is small when q^ is small. Our next step involves the 
specification or computation of Iraf-,^{s) by exploring the fact that only the pn state which couples to 
the pion contributes to Imf-,^{s) in the momentum transfer region of interest to this paper. An one loop 
computation of Imf-,^{s) is known to be insufficient to account for experimental data in other studies 



21, wM, the pvr correlation, which forms a resonance near 1.3 GeV, is required. We therefore propose 



tEe following form for /m/^(s) 



ImUq') = l^'^-p.,Ae) (4.5.23) 

'ignNN 47r 



with the reduced density of state 



P^Mq^) = Po{Q^)i^ + T-2 \l^l-2f 2U2 1 (4.5.24) 



and 



-(2-, /I {mp + m^Y {rrip-m^f ( ml -ml 
= q^ V V2 V~^^ 



-2 {mp + m^)\ (4.5.25) 

where 6{x) is the step function with a value of unity for positive x, sr is chosen to be 1.69 GeV^, A 
characterizes the strength of the /jvr resonance in Imf.,^{q^) and Pq{sr)Ib characterizes the width of the 
resonance. 



The above form is chosen so that when A = 0, ImfT^{q^) is the one loop result in the Feynman-t' 
Hooft gauge (for the p propagator). The /otttt interaction piece of Lagrangian density used for evaluation 
of Imf-,^[q^) is 

Cp^ = 5p..e"''=^"5''^V^- (4.5.26) 

^The orthogonality relation between vertex functions in the relativistic case may need generalization. It is however 
expected to exist and to be unique. So the projection procedure is not ambiguous. In many analytic diagramatical 
calculations based on simple models, the pion contribution term can simpUy be read out from the three point function. 
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The corresponding piece of the axial vector current operator, which can be obtained from the Noether 
theorem, can be written as 

= 5p..^6'^'^vrV^, (4.5.27) 

QtvNN 

where use has been made of the hnear a model relation rriN = QirNNcr- 

The value for A, AttAtat and Ib is adjusted so that the minimum value of the following function 

f{X,A^,.,lB) - ^^^^^^ e^., (4.5.28) 

with q1 = q^in + ^qmax ~ qmin/^ ^ achieved. The value of is chosen to be 100. = —0.6 GeV 
and q^ax ~ ^-2 GeV. The factor e^"? is used to put more weight on small \q^\ region where the fit tends 
to be poor. 

In all cases studied and presented in table 4J, A^attv < 0.95 GeV. If a value g'^^/Air = 1.0 is taken. 



Table 4.1: The results of fitting, where c/a = 1-26, g^NN = 13.4, U = 93.2 MeV, g^j„ = -0.6 GeV^ and 
ql„^ = 0.2 GeV^. The unit for all but A and gpniT is GeV. A and gpmr are dimensionless. The quantities 
with a star on top is chosen by physical considerations. 





^ttNN 


A 


VIb 




Ml 


V^th 


1.0 


0.94 


1.12 


1.12 


1.30 


1.00 


0.91 


1.5 


0.93 


0.41 


1.26 


1.30 


1.00 


0.91 


2.0 


0.91 


-0.0023 


1.34 


1.30 


1.00 


0.91 


2.5 


0.89 


-0.25 


1.27 


1.30 


1.00 


0.91 


2.9 


0.88 


-0.39 


1.30 


1.30 


1.00 


0.91 


3.5 


0.86 


-0.56 


1.35 


1.30 


1.00 


0.91 



the qualitative shape of ImfT^{(f') in Fig. 4T obtained by a minimization of Eq. 4.5.28 is similar to the 
ImT{(p') of Ref. [^2|. Quantitatively, it has a broader width. The phenomenological value for ^ptttt can 
be deduced from the /> — > vrvr decay process. It has a value satisfies fl'pTrTr/^Tr ~ 2.9. Using this value of 
gpT^^, ImfT^{q^) is obtained by minimization. The result is given in Fig. [4.1[ It's drastically different 
in shape from that of ImT{q^) in Ref. |22]. In fact, instead of increasing the density of states relative 
to t he one loop result, the resonance contribution decreases the density of states in order to satisfy 
Eq. 4.3. 15| . The reduction of density of state indicates either the solution is unphysical (for a normal 



resonance) or there is a competing resonance in another channel that couples to the pionic channel we 
are dealing with. But what can the "other resonance" channel be? There is no known resonance there. 
There seems to be an inconsistency. 

Of course this conclusion has to be checked by other computations of /?Ti/,r(^^) in different model 
Lagrangians or in more fundamental ones like the lattice QCD calculation. It is a worthy topic to be 
examined in the future. 
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Figure 4.1: The pvr reduced density of states in the pionic channel as a function of obtained from 
the fitting procedure for two values of the pvrvr couphng constant. For comparison, the one loop (A = 0) 
reduced density of state is drawn with a dotted line. Here o-p-K-K = Qp-jx-K 

/47r. 



4.6 Extended PCAC relation for a nucleon 



Two relationships between experimentahy accessible nucleon form factors gA{(f), 9p{(f) and g-wNNiff') 
and the less well known fnig'^) off the pion mass shell are established based on the facts that 1) the 
QCD Lagrangian has an approximate S U(2)l X SU (2)/j symmetry, which is explicitly broken down only 
by a small current mass term (see Eq. 4.1.1| ) 2) this chiral symmetry is spontaneously broken down 
to an flavor SU{2)v symmetry with pion, which is lighter than it should be as an ordinary hadron, as 
the Goldstone boson of the symmetry breaking 3) the pion, being lighter than other normal hadronic 
particles, should dominate the low mome ntum transfer reactions in certain channel, which is expressed 
as the PCAC relation given by Eq. 4.1.4. 



The first relation Eq. 4.3.14 passed the experimental test using the available data. The second one 
seems to be in trouble when it is confronted with our empirical knowledge. In general, if the deviation 
is finally established, the remedy for it can be either at a fundamental level, which means to modify 
the contemporary field theoretical framework or QCD, or at the structural level in the sense of revising 
our notion of what the structure of the physical system under investigation is. We adopt the second 
alternative in the following since even if the unlikely possibility that the inconsistency discussed in the 
above section is finally elliminated after more detailed study, it is still an intellectual challenge how 
the PCAC relation for a nucleon can be extended within the current theoretical frame work. After all, 
it is hard to imagine the portion of the vacuum state inside a nucleon remains unchanged under the 
influence of the baryon density of the valence quarks that are compacted into a region of order 1 fm. 

The basic idea of the c urrent approach to extend the PCAC relation consists of a modification of the 
correspondance Eq. 4.1.4 when its matrix elements are taken between nucleon states by assuming that 
besides the pion, there is a different set of soft modes inside a nucleon originating from the spontaneous 
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breaking of the chiral SU{2)l x SU{2)ji symmetry down to the same flavor SU{2)v symmetry. Such 
a possibihty is theoretically investigated in Refs. ^ and are presented in the first p art. I f there are 



such a set of these soft modes inside a nucleon, then the left hand side (l.h.s.) of Eq. 4.1.4 in the low 
momentum transfer regions is saturated not only by the pions, but also these soft modes (Goldstone 
diquark excitations), which exist if a localized /3-phase is assumed. 

These color carrying soft modes are confined inside the nucleon and therefore c an not be directly 



observed like the pions but can has an effect by have non-pole contributions to Eq. 4.3.15 . The result 




IS 

„2 



mgAiq') = 9.NN{q')U{<l') + {<1- Ovill, (4.6.29) 

where ??(g^) is a function related to the product of the strength of the Goldstone diquark excitations 
and their propagator, which has no pole at low energies (see Ref. ||8| for a more careful discussion on 
this point) due to the fact that they are confined inside the nucleon. It is worth mentioning that the 
additional term has no effects on the pion mass shell — m\ due to the — factor. 

With this term added, the potential inconsistence between theory and data discussed in the previous 
section can in principle be resolved. It is clearly interesting to further investigate the existence or the 
extend of the inconsistency discussed using various other models and computation procedures. 
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Chapter 5 



Gerasimov— Drell— Hearn Sum Rule for a 
Nucleon 



The Gerasimov-Drell-Hearn (GDH) sum rule 26| is a relation between the difference of two polarized 



total photon- spin-1/2 Dirac particle cross sections and the corresponding particle's anomalous magnetic 
moment. It can be expressed as 



m 



where ccem = 1/137 is the fine structure constant, (T3/2 is the total cross section when the photon helicity 
is in the same direction as the target particle's spin polarization, ai/2 is the one when the photon helicity 
antiparallels to the particle's spin polarization and k is the anomalous magnetic moment of the particle. 

It can be derived if the corresponding forward photon-particle scattering amplitude decreases fast 
enough at high energies so that an unsubtracted dispersion relation for the amplitude can be written 
down. For the case of photon-nucleon scattering, the large energy behavior can be inferred from 
phenomenology, which tells us that the GDH sum rule ought to be satisfied. The current confrontation 
of the GDH sum rule for a nucleon to known experimental data from pion photo-production result in 
disagreement p^, pi]. What causes the discrepancy is not yet clear. 



Maybe the experimental information used to saturate the l.h.s. of Eq. p.0.l| is not enough. It 
may also be that there is a true need of an extension of the GDH sum rule on the theoretical part. 
This study focus its attention to a theoretical study of the possibility of an extension of the GDH sum 
rule consistent with Lorentz covariance and local electromagnetic (EM) gauge invariance using the field 
theoretical methods. We shall review briefly various issues concerning the modiflcation of the GDH sum 
rule proposed. Details are given in Refs. B 



5.1 Photon— nucleon compton scattering 



The photon-nucleon forward Compton scattering amplitude is related to the following covariant one 
photon irreducible correlator between EM current operator J'^(x) 

T^'^p,q) = i j (fxe''i-''{pS\T*J>'{x)r{Q)\pS), (5.1.2) 
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where \pS) is a nucleon state with 4-momentum pf^, polarization S^^ and T* represents time ordering 
with proper Schwinger terms added. Taking into account of the conservation of parity in the EM 
interaction, the amphtude can be parameterized by eight invariant amphtudes -F1...8 the following way 

+ iF^ {fa'^'^q^ - p^'a^'^q^) + iF, {q^a'^'^q^ - q^'a^'^q^) + iFge^^^^^p^^T'] U(j>S). 

(5.1.3) 

The invariance amplitudes -F1...8 are functions of and v = p ■ q/m with m the mass of a nucleon. 

Since the amplitude T^" satisfies the Ward identity 

qt.T''''{p,q) = (5.1.4) 

due to the gauge invariance and commutativity of the EM charge density operators at equal-time, it is 
usually written in a reduced form, namely. 



mu ,,\ ,, mv 



T^^p^q) = S,[-g^^+^-^^+S,[p^-^q^][p'^-—q^ 

-l2 



with 



Mie^-"^„5^ - imvA^e^^'^f'q^ [Sp - ^Pp) (5.1.5) 

Si{q\y) = -Fiiq\u) = -q^F2iq\u)-miyF4{q^u), (5.1.6) 

^2(g',^) = Fs{q^i^) = ^F,{q^u), (5.1.7) 

mu 

Ai{q\u) = mFe{q\u) + uFs{q^u), (5.1.8) 

A2iq\i^) = -F^{q^,v) + F^{q^,u). (5.1.9) 
m 

The forward cross section of the compton scattering is related to the amplitude Mfi given in to 
following 

Mfi ~ e'^e.T'^^ (5.1.10) 

where e' and e are the final and initial photon polarizations. In the nucleon rest frame, the polarization 
dependent part of the forward cross section depends only on Ai when the nucleon polarization direction 
= (0, S) is along direction of the photon propagation. 



The invariant amplitude Ai enjoys a sum rule given by Eq. |5.0.1| . Its derivation involves three issues 
discussed in the following 1) is the use of the infinite momentum frame legitimate? 2) can the high 
energy {v) behavior of Ai{q^ , v) be estimated using the Regge asymptotics? 3) what's the role of gauge 
invariance in determining the large v behavior of Aii^q^, v)l 



5.2 Sum rules and infinite momentum frame 

The need of using an infinite momentum frame in deriving a fixed q^ sum rule like the GDH sum is 



discussed and emphasised in Ref. [29|. It's validity in the derivation is generally assumed. It is also 
assumed in this work. 
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5.3 Regge behaviors and physical states 



The large energy physical hadronic scattering amplitudes follow certain power law (in energy) according 
to Regge asymptotics iQ, namely 

lim Ar^u" (5.3.11) 

where v is the energy variable and a is the leading Regge trajectory of the corresponding channel. The 
highest trajectory is that of the Permeron with an q(0) = 1 [p^]. Since the Permeron has a spin one (in 
the forward direction), it can not be exchanged in some channels of reaction due to helicity conservation; 
this allows the corresponding amplitude to change slower than v. If the change is slow enough, it allows 
certain super convergence relation to be derived for it in the corresponding channel for physical hadronic 
scattering amplitude 



The problems of using Regge asymptotic behavior in deriving sum rules for matrix elements of a 
current (between physical hadronic states) is known from current algebra studies |32, 33 1, in which 



it was found that a hypothetical case of isovector photon-nucleon scattering there appears the need 
of a "J=l fixed pole" in the amplitude for the matrix elements of the time ordered current-current 
correlator. This problem can be understood if one realizes that the matrix elements of the time ordered 
current-current correlator does not correspond to a 4-points physical hadronic scattering amplitude 
directly. The phenomenological Regge asymptotics may fail for these amplitudes. This situation opens 
the door for a modification of those sum rule that are derived under the assumption of Regge behavior 
for the matrix elements of the time ordered current-current correlator. For the GDH sum rule, this is 
suggested in Refs. [p^, |35[. 

For the realistic photon-nucleon scattering, this kind of "J=l fixed pole" effects can not be straight 
forwardly introduced in such a way that the results are consistent with current algebra. This is because 
it causes a modification of the commutation relation between the time component of the EM current 
operator at equal-time, which tells us that they commute with each other. This requirement, together 



with gauge invariance, are expressed simply by Eq. |5.1.4| . A modification of the commutation relation 
between the charge density operators at equal-time is suggested in Refs. |3f:, 5^]. It leads to a modifica- 



tion of the Ward identity Eq. |5.1.4| , so that a mutual consistency of the arguments can be maintained. 
Whether or not such a way of modifying GDH sum rule actually is consistent with phenomenology 
remains to be shown however. 



5.4 Gauge invariance constraints 



The scheme proposed in this study is different in that the GDH sum rule is modified by introducing 
an effective "J=l fixed pole" effect without modifying the commutation relation between the charge 
density operators at equal-time. This is achieved by assuming there is a localized breaking down of the 
U{l)em gauge symmetry inside a nucleon The "J=l fixed pole" effects are provided by the would-be 
Goldstone boson of the symmetry breaking, which does not belong to physical excitation spectrum for 
a local gauge symmetry p8|. 



From the Ward identity Eq. ^.1.4| , there are three constraints for -F1...8 in Eq. ^.1.3 , namely, 

Fi{q'',i^)-q^F2iq^u)-muF4{q'',i^) = 0, (5.4.12) 
miyF3{q^u)-q^F4{q^u) = 0, (5.4.13) 
F5{q\u)-muFe{q^u)-q^Fr{q\u) = 0. (5.4.14) 
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Eq. 5.4.14 is relevant to the possibility of modifying the GDH sum rule. 



The large v behavior of F5{q , v) can be determined to be [H] -F5 ~ . If the following assumption 



is made, namely, 



lim q^Fj = 0, 



(5.4.15) 



the n the la rge v behavior of at = is following Eq. 5.4.14. Under the assumption given by 
Eq. 5.4.15, the large v behavior of Ai, which is relevant to the GDH sum rule, is controlled by that of 
-Fg- The large v behavior of F^ is bounded by the Regge asymptotics because it is an gauge invariant 
invariant amplitude by itself, which means that it connects to physical scattering amplitude only (for 
a detailed discussion on this point, see Refs. [Q, |3^). Thus Fg ~ ~ with a < 1 from helicity 
amplitude analysis |4^. Thus gauge invariance requirement expressed by Eq. 5.1.4 together with Eq. 
5.4.15 eliminates the possibility of a modification of the GDH sum rule. 



5.5 The possibility of modifying GDH sum rule 



The above analysis show that if GDH sum rule is to be modified in a way that respect the local EM gauge 
inva riance and the Regge asymptotic behavior for physical hadronic scattering amplitudes, assumption 
Eq. 5.4.15 has to be relinquished by letting F-j to has a pole like behaviorQin at large u. This implies 



that the U{l)em symmetry corresponding to EM is spontaneously broken down inside a nucleon 
More specifically, introducing an order parameter p^o, defined as 

q^Fj 

Poo = — lim lim (5.5.16) 

for the EM gauge symmetry inside a nucleonpl. If the U{l)em gauge symmetry is spontaneously broken 
down, then there is a massless pole in Fj in tne u ^ 00 limit so that poo can be nonzero. With pao, the 
GDH sum rule is modified to 



f 

Jo 



du = fL^(^n'' + 2m^Poo). (5.5.17) 



Numerically, the value for and p^ and their difference can be extracted from the results of 



the integration of the total photon-nucleon cross sections given in Ref. |28]. They take the value 
p^ = 2.9 X 10~^/m^ and p^ = —2.5 x 10~^ fm?. As a final remark for this section, it should be 
mentioned that p^" has a dimension of length squared. It can be written as pg^" = 1/A^. It is 
interesting to note that A ~ IGeV, which correspond to the natural scale of the spontaneous chiral 
symmetry breaking in strong interaction. 

It shows that if the current discrepancy between cross section (T3/2 and (T1/2 obtained from the the 
photo pion production on a nucleon data and the GDH sum rule is genuine, then a localized spontaneous 
breaking down of the U{l)em symmetry inside a nucleon has to be introduced. As it is discussed in part 
I of this report, such a symmetry breaking is in principle possible since the U{l)em gauge symmetry is 
spontaneously broken down in the /3- and uj- phases |5|, |^ . 



^Ref. [g] provides a more precise meaning for this statement. 

^The first fimit — > 00 is necessary since locaUzed spontaneous gauge symmetry breaking is discussed here. A finite 
region in space only looks more and more like an infinite system when smaller and smaller distances or higher and higher 
energies are probed. 
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Chapter 6 



Rapidity Correlation of BB in High 
Energy e^e~ Annihilation 



The hadronization processes in the e^e^ annihilation can be described reasonably well by a chain like 
picture in terms of string (or flux tube) fragmentation. For the meson production, quark and antiquark 
pair is created by a breaking of the string that connects the parent quark-antiquark pairs. At the 
breaking point, any nonvanishing conserved quantum number like charge, baryon number, etc. are 
created in conjugate pairs so that they cancel locally (in space-time) there. The reason behind it is that 
the string is considered to be neutral in these quantum numbers. Monte Carlo simulation programs like 
the JETSET [^] and Herwig |42|, which are based upon such a picture describe the experimental data 
well. 

This picture of hadronization implies strong rapidity correlation between pairs of hadrons that 
conjugate to each other in the multihadron final state. Since, before the string fragmentation, it is in 
a stage of uniform expansion, the neighboring hadrons produced on the string by a fragmentation has 
rapidity close to each other than others. It is difficult to test the consequences of this picture using the 
mesonic component of the final state since 1) the correlating meson pairs are hard to identify and 2) the 
interaction between the mesonic particles are stronger than that of the baryons in the e^e~ annihilation. 
The interaction distorts the rapidity information of the meson at the time when it is created. On the 
other hand, larger than 50% of the baryons escape the hadronic clouds before interaction effects grow; 
they therefore carry the rapidity information at their production time One of the best characters 
of the string fragmentation picture for hadronization in the e^e~ annihilation processes is the strong 
rapidity difference correlation between a pair of baryon and antibaryon {BB). 

The experimental examination of it was carried out by observing the AA rapidity difference cor- 



relation |10| in e+e" annihilation at y/s = 91 GeV by OPAL collaboration. A "surprise" was found 
that the width of the rapidity difference correlation between AA pairs are wider than it is expected 
from the one predicted by the above mentioned fragmentation model, which consists of the creation of 
diquark-antidiquark pairs on the string for the formation of baryons, similar to the quark-antiquark 
pair creation for mesons production. In order to describe the data, the so-called popcorn mechanism 
has to be introduced, and, in addition, the percentage of the popcorn configuration has to be large 
in the fragmentation processes |l^ . 

The popcorn mechnism implies a non-local production of conjugate baryon number pairs on the 
string under the fragmentation since, in between them, there has to be a finite space that contains the 
meson. The conclusion that conjugate baryon numbers carried by diquark and antidiquark on the string 
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are not created at the same spacetime point is hard to understand if the string that fragments into BB 
pairs are neutral in baryon number. According to the principle of relativisty, it violates the classical 
causality which requires local cancellation of baryon number and antibaryon number at the time when 
they are created from a baryon neutral source; the speed of their separation should be less than the 
speed of light. 



This difficulty can however be solved if one assumes that there exsits /3- or uj- phase inside the string 
under the fragmentation. The reason, as it is discussed in section 3A of part I is that, spontaneous 
separation of baryon numbers are favored in these phases (see, for example, Fig. 3.3). Therefore if one 
of the above mentioned phases is present inside the string before a fragmentation, the baryon density 
on the string can be non-zero at a specific point; it is positive or negative alternatively along the string. 
Thus, by assuming the existence of the /?- or uj- phase discussed in part I of the report, we can gain a 
deeper understanding of the empirical popcorn mechanism needed in phenomenology. 



Before closing this short chapter, let's reminded the reader that this finding, if proved genuine 
in the future after more comprehensive studies, could provide a domestic experimental basi^ for the 
baryogensis mechanism in a matter-antimatter symmetric universe [^] under the standard cosmology 
(big-bang), which is based upon the spontaneous separation of baryonic matter and antibaryonic matter 
at certain time during the evolution of the univers [|, without violating the observational constraints 
p^. It is certainly a worthy topic to be further studied. 



In the sense that the physical processes underlying the mechnism can be created and studied here, on earth. 
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Chapter 7 

Summary and Outlook 



The theoretically possible phases in an interacting massless two flavor quark system are discussed in the 
chapters of part I of this report by introducing model Lagrangians. Four phases are found. Some of their 
properties are discussed while others are mentioned with references containing more details provided. 
In part II, three different observations, which is currently considered difficult to understand using the 
conventional picture, are discussed in the light of the findings presented in part I of this report. We 
find it likely that localized diquark condensed phases exist inside a nucleon and hadronic excited states 
(within the flux tube). 

In the heavy ion collisions, baryon number density and the volume containing it can be much larger 
than that of a nucleon. It is thus very probable that the (3- or lo- phase can be produced despite the 
short time duration of the collision. What signals their existence existence and how to find them are 
questions that can be further investigated. 
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